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ABSTRACT 

Using results from radio and X-ray observations of millisecond pulsars in 47 Tucanae, and extensive 
HST U, V, I imaging of the globular cluster core, we have derived a common astrometric solution good to 
< O'.'l. A close positional coincidence is found for 47 Tuc U, a 4.3ms pulsar in a 0.429 day orbit, detected 
in radio and X-rays, with an my = 20.9 blue star. Analysis of extensive time series data for this optical 
candidate shows a 0.004 magnitude semi-amplitude variation at the period and phase expected from the 
radio ephemeris, and the optical variations are spatially coincident with the candidate. This provides 
secure optical detection of the white dwarf companion to the millisecond pulsar, the first such detection 
in a globular cluster, allowing for comparisons to recent models for such companions with dependencies 
on mass and age. 



Subject headings: binaries: general - globular clusters: individual (47 Tucanae) 
(PSR J0024-7203U) - pulsars: general 



pulsars: individual 



1. INTRODUCTION 

The huge improvement in spatial resolution of Chan- 
dra over previous X-ray missions is proving, like Hubble 
Space Telescope (HST), to have significant impact on glob- 
ular cluster astrophysics. For example, using Chandra, 108 
sources have been detected in the central 2' x 2'.5 of the 
massive globular cluster 47 Tucanae (NGC104; Grindlay et 
al. 2001), compared to the detection of 9 central sources 
by ROSAT (Verbunt & Hasinger 1998). These Chandra 
sources are a mixture of millisecond pulsars (MSPs), cat- 
aclysmic variables (CVs), quiescent low-mass X-ray bina- 
ries, and active binaries, many of which have also been 
detected (Edmonds et al. 2001, in preparation) in our ex- 
tensive 8.3 day _ffS'T/WFPC2 observations of this cluster 
(Gilliland et al. 2000). Fifteen of the 20 MSPs detected in 
47 Tuc by Camilo et al. (2000) with the Parkes radio tele- 
scope have timing positions accurate to <0'.'01 (Freire et 
al. 2001; hereafter FCLOl), and 10 of these are detected in 
X-rays with the CIAO program wavdetect (errors < 0'.'2), 
allowing for excellent astrometry between the radio and X- 
ray coordinate frames. Combined with the large numbers 
of CVs and active binaries detected with both Chandra 
and HST, we exploit here the rare opportunity to perform 
astrometry at the ^^0.1" level between X-ray, optical, and 
radio observations, using only relative astrometry. 

Of the 8 binary MSPs in 47 Tuc with timing solu- 
tions (FCLOl), 5 have companions with masses of ^0.15- 
O.2OM0 (Camilo et al. 2000); by analogy with field MSPs 
these are expected to be helium white dwarfs (He WDs). 
Using our precise astrometry we report here the unam- 



biguous detection of a He WD companion to one of these 
MSPs (PSR J0024-7203U, hereafter 47 Tuc U), the first 
such detection in a globular cluster. (He WDs have been 
reported in the globular cluster NGC 6397 by Cool et al. 
1998, Edmonds et al. 1999, and Taylor et al. 2001, but 
none of these have yet been associated with an MSP.) The 
astrometry, photometry and time series for this identifi- 
cation are described here, along with a discussion of their 
consequences for the age of 47 Tuc U, the nature of He 
WDs and their cooling, and the heating effects of an MSP 
in a close binary. 

2. OBSERVATIONS AND ANALYSIS 



2.1. Astrometry 

We assumed the astrometric coordinate system of the 
MSPs (based on the JPL DE200 planetary ephemeris; 
FCLOl) as our reference. First, we transferred the Chan- 
dra coordinate frame onto the radio frame using only the 
relatively isolated MSPs with accurate Chandra positions 
(six MSPs including 47 Tuc U). We calculated and sub- 
tracted the mean positional offsets between the X-ray and 
radio coordinate systems and found the rms residuals were 
0'.'136 and 0'.'037 for RA and Dec respectively (after cor- 
recting for small linear correlations between RA and RA 
error and Dec and Dec error). The total error of 0'.'14 
compares favorably with the mean predicted wavdetect 
error, at the relevant count levels and in the absence of 
background and crowding, of 0'.'15. 

We then fixed the HST coordinate system, with posi- 
tions based on the STSDAS program metric, onto that of 
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Chandra using the X-ray detection of 6 CVs and active bi- 
naries found on the WF2 chip of the Gilhland et al. (2000) 
HST/WFFC2 data (this chip contains 4 of the 6 MSPs, 
including 47 Tuc U, used to perform the radio to Chandra 
astrometry). The rms residuals were found to be 0'.'028 
and 0'.'079, for RA and Dec respectively. Then, assuming 
the errors are random and dividing by \/6, the resulting 
la errors matching X-ray to optical are O'.'Oll and 0'.'032 
for RA and Dec, and for radio to optical are 0'.'057 and 
0'.'049. 

Using this astrometry we searched directly for MSP 
companions in the HST images and found an excellent 
candidate for a companion to 47 Tuc U (hereafter Uopt) 
whose position is given in Table 1. Finding charts for Uopt 
(Fig. 1) show the [/-band and F-band images, plus the 
3(7 error circle for Wll (the Chandra X-ray source corre- 
sponding to 47 Tuc U) and the 3cr error ellipse for 47 Tuc 
U. The offset between Wll and Uopt is 0'.'055, and between 
47 Tuc U and Uopt is O'.'IO. 

2.2. Photometry 

Figure 1 shows that Uopt is clearly blue, since other stars 
in the V image of comparable brightness to Uopt are barely 
visible in the U image. This impression is confirmed by 
our U vs U — V and V ys V — I color-magnitude diagrams 
(CMDs) from PSF-fitting photometry shown in Figure 2 
(we have transformed from iJST filters to Johnson-Cousins 
using Holtzman et al. 1995; see also Table 1). The star 
Uopt (labeled 'U'; My = 7.47) is blueward of the main 
sequence (MS) in both CMDs, where the internal error for 
Uopt is < 0.03 mag for each filter, with absolute errors of 
~ 0.05 mag in V and / and ~ 0.1 mag in U, including 
transformation uncertainties. A likely CV is labeled ('1') 
as are 5 blue stars (circles) confirmed by visual examina- 
tion to be relatively bright in U like Uopt- One of these 
blue stars is probably a CO WD (see Zoccah et al. 2001), 
two others are variable and likely CO WD-MS binaries 
and the other two are possibly also binaries. Visual exam- 
ination showed that the other apparently blue stars, with 
V < 22, are caused by crowding and photometric error, 
particularly in the crowded portions of the image near- 
est to the cluster center. Many of the fainter stars with 
U > 22 and < C/ — U < 1 are bona- fide blue stars, likely 
to be dominated by MS stars from the SMC (see Zoccali 
et al. 2001). 

The very small number of stars that are blue and de- 
tectable in each of U, V and / highlights the unusual na- 
ture of Uopt and its likely association with 47 Tuc U. With 
7 such blue stars the probability (averaged over the whole 
WF2 image) of detecting one in a circular aperture of 0'.'15 
radius is only 9x10^^. 

2.3. Time Series 

The HST/WFPC2 data of Gilhland et al. (2000) contain 
636 U-band and 653 /-band 160 s images evenly spread 
over 8.3 days. Typical rms noise values in fractional in- 
tensity (Gilhland et al. 2000) are - 0.003 at U = 17 and 
^ 0.013 at U = 20. Since Uopt is relatively well isolated, 
the quality of its V and / time series, based on its ex- 
tracted aperture of 0'.'34, is excellent. Its rms noise values 
are 0.021 (U) and 0.034 (/), or -0.1 and 0.6cr away from 
the mean rms values for the 666 (V) or 551 (/) stars within 



0.25 mag of Uopt, after iteratively removing > 3a outliers. 

We then searched for evidence of low amplitude varia- 
tions in Uopt- Lomb-Scargle (Scargle 1982) power spectra 
of the V and / time series for 47 Tuc U are shown in the 
upper panel of Figure 3. This shows a peak in the U-band, 
at (2.681±0.07) x lO^^Hz, which is coincident with the or- 
bital period for 47 Tuc U (2.697 x lO^^Hz) determined 
from radio observations. Although this peak is not statis- 
tically significant in a blind search, only one other peak, 
when measured out to 0.45 niHz (the 'average' Nyquist 
frequency; Scargle 1982) has a higher signal. No corre- 
sponding signal is seen in the /-band. 

We then carried out a least-squares fit of a sinusoid to 
the time series, by coordinating the radio ephemeris with 
the optical and correcting the mid-exposure HST times to 
the barycentcr, using the IRAF task rvcorrect. We de- 
fine phases [T - (Tasc - /'/4)]/P (where T is the barycen- 
tric time, Tasc is the time of ascending node, when the pul- 
sar is moving with maximum velocity away from us, and P 
is its orbital period) so that when phase = 0.0 the heated 
hemisphere of Uopt is facing towards us, giving maximum 
light (in a model where Uopt is heated by the pulsar wind) . 
The derived period, time of optical maximum (phase = 0) 
and semi-amplitude for the U-band time series are given 
in Table 1 {V phase plots are shown in the lower-left panel 
of Fig. 3). For the 'unforced' solution the period is 0.7a 
away from the orbital radio period, and the phase is 1.9cr 
away from phase zero for the radio observations (as defined 
above). The amplitude is significant at the SAa level (the 
significance of this measurement is > 3Aa when the ~1/10 
chance probability of the radio and optical phases align- 
ing is considered) . By forcing the period and phase of the 
sinusoidal variation to equal the radio values (the 'forced' 
solution), an amplitude of 0.0033±0.0012 is derived, sig- 
nificant at the 2.8(7 level; a forced fit in the /-band yielded 
0.0013±0.0018. Marginal evidence is found for variation 
in the much more limited (26 images of length 160 s - 900 
s) C/-band data (see Fig. 3). 

To test whether this small variation in V is intrinsic to 
Uopt rather than a larger variation from a nearby fainter 
star we examined the power spectra of 3 fainter stars (with 
V > 22) lying within ~ 0'.'5 of Uopt- None showed signif- 
icant variations. We also examined a 'phase difference' 
image derived by adding together images where the phase 
of Uopt is within 0.07 of 1.0 and subtracting similar im- 
ages for phase values of 0.5±0.07. This procedure gave a 
stellar-shaped source coincident, at the 0.5 pixel (0'.'05) er- 
ror level, with Uopt , exactly as expected if Uopt is the only 
star varying with the radio period and phase. We there- 
fore conclude we have made a solid detection of variability 
from Uopt at the expected period and phase, confirming it 
as the 47 Tuc U companion. 

3. DISCUSSION 

Neutron stars (NSs) in binary radio pulsar systems have 
average masses of Mns = 1-35 ± 0.05 M© (Thorsett & 
Chakrabarty 1999). For 47 Tuc U, considering a NS mass 
1.30 M0 and an edge-on orbit (inclination i ~ 90°), we de- 
rive from the mass function a minimum mass for Uopt of 
Mu — 0.120 Mq. This 4.3ms pulsar has an approximate 
characteristic age of t ^ 2 Gyr and spin-down luminosity 
of iJ ~ 2 X 10"^"* erg s^^ (Grindlay et al. 2001, in prepara- 
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tion) . Such recycled pulsars may have accreted significant 
amounts of mass, and for every additional 0.05 M© for the 
NS, Mu's minimum value increases by 0.003 M0. 

We now compare the U, V, I colors of Uopt with expec- 
tations from H and He WD models to test whether it 
is a plausible He WD and if so, to estimate TeS, log g, 
and the mass and age of this star. We focus on compar- 
isons with the models of Serenelli et al. (2001; hereafter 
SAROl), who conveniently provide broadband color infor- 
mation. They find that for masses ranging from 0.18 M0 
to 0.41 Mq, diffusion-induced hydrogen-shell flashes take 
place which deplete the H envelope and cause evolution to 
proceed relatively quickly, while the He WDs with mass 

< O.I8M0 evolve much more slowly because they retain 
relatively thick hydrogen envelopes. Althaus, Serenelli, & 
Benvenuto (2001) and SAROl show that this dichotomy 
between low and high mass He WDs removes the discrep- 
ancies between MSP characteristic ages and WD cooling 
ages for several field systems including PSRs J0034— 0534, 
J1012-h5307, J1713-h0747 and B1855-h09. 

Figure 2 shows the He WD cooling curves of SAROl for 
masses of 0.169 - 0.406 Mq. Since Uopt hes well to the 
red of the WD cooling sequence of Bergeron, Wesemael, 
& Beauchamp (1995) in both CMDs, it is unlikely to be 
a ~0.55 M© CO WD. The lowest mass models of SAROl 
offer much better agreement, and since Uopt falls slightly 
to the red of the lowest mass sequence we infer that My 

< 0.17Mq (implying i > 45°). This would mean that 47 
Tuc U is similar to the field system PSR J1012-h5307, with 
period = 0.605 d and mass ~ O.I6M0 (Callanan, Gar- 
navich, & Koester 1998). 

The 0.1 7M0 SAROl model with the same absolute mag- 
nitude as Uopt has Toff — 11000 K, log g — 5.6, luminosity 
= 0.14 Lq, a radius^ of ~ 7.3 x lO^cm and an age of 
^ 0.6 Gyr. This cooling age is somewhat less than the 
characteristic age estimated by Grindlay et al. (2001, in 
preparation) and is at the low end of the conservative lim- 
its set by FCLOl of 0.4 < r < 4.2 Gyr. Only the lowest 
mass model of SAROl remains relatively bright at these 
significant cooling ages, while all of the higher mass WDs 
reach much fainter levels (My ^ 12) at a cooling age of 
~ 1 Gyr (see Fig. 2). Our data therefore provides tenta- 
tive support for the evolutionary calculations of SAROl, 
in particular for the dichotomy in ages between different 
mass He WDs. For comparison, the He WD models of 
Hansen & Phinney (1998) (for all masses) reach M/=11.5 
or fainter after 1 Gyr, for a relatively thick hydrogen en- 
velope of 3 x 10"" Mq. 

Accretion should have ceased in the 47 Tuc U sys- 
tem, to explain the faint X-ray luminosity (2 x lO'^^crg 
s^^; Grindlay et al. 2001, in preparation) and the 
constant, uncclipsed radio emission (FCLOl). Using 



the Roche-lobe formula from Paczyhski (1971) {r/a = 
0.462[(Mu/(Mns + Mu)]^/^ where r is the Roche-lobe ra- 
dius, and a is the binary separation), the WD radius given 
above (for Mu = 0.17 Mq; i = 45°) and the binary sepa- 
ration from Kepler's Third Law, Uopt should under -fill its 
Roche lobe by a factor of ^6. Therefore, as expected, no 
accretion will be occurring. 

Assuming that the MSP radiates its wind isotropically, 
and that i = 45°, then the E estimate implies that the 
pulsar energy intercepted by Uopt should be ~ 7 x lO^^erg 
s~^ or O.OlSx the luminosity of Uopt (assuming that the 
spin-down energy is re-radiated as a black-body, we derive 
a temperature of ^ 5200K, peaked at ^ 5600A, for the 
re-radiated emission). Therefore, assuming that only one 
side of Uopt is heated we should see a relative variability 
amplitude of ^ 0.005, allowing a factor of 1.4 decrease in 
possible amplitude because of the inclination. This is only 
slightly larger than the amplitude of the detected V-band 
variation, implying a high efficiency factor for re-radiation 
of the pulsar's wind as optical emission. However, because 
the variability expected from the effects of heating is small, 
we infer that the possible (2.3(7) much larger amplitude 
variation in the C/-band is either spurious or evidence for 
a separate source of emission. 

4. SUMMARY AND FOLLOWUP OBSERVATIONS 

We have identified the binary pulsar 47 Tuc U with an 
optically variable blue star with period and phase of maxi- 
mum light in excellent agreement with that expected from 
the precise radio ephemeris (FCLOl). The companion is 
likely to be a He white dwarf with My ^0.17 Mq, based 
on the published models of SAROl. 

Although the optical time series photometry reported 
here is unlikely to be improved upon soon, time resolved 
spectroscopy with _ff5'T'/STIS of Uopt appears promising. 
The co-added Doppler corrected spectrum could be used 
to spectroscopically determine Toff and log g, and hence 
Mu (see e.g. Edmonds et al. 1999). Combining My with 
the radial velocity amplitude of Uopt (Ku) and the mea- 
sured Kns from the radio measurements would give a di- 
rect measurement of Mns • 
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Table 1 
Positional, photometric and time series information for 47 Tuc U and its optical companion Uopt 



(J2000) 


Dec'^ 
(J2000) 


U 


V 


/ 


period(y) 
(days) 


T(phase=0)'' 
(MJD) 


amp(y) 
(mag) 


00 24 09.8325(5) 


-72 03 59.667(3) 


20.6(1) 


20.91(5) 


20.70(5) 


0.432(4) 


51367.32(2) 


0.004(1) 



'^Coordinates from Freire et al. (2001) 
''Time of optical maximum (see text) 
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^^ 4^^ ^ U image 




Fig. 1. — Finding chart for 47 Tuc U showing the co-added (over 8.3 d), oversampled images for the U- and V-band. The solid hne shows 
the 3cr error circle (radius = 3 X 0''14) for the X-ray source Wll, and the dotted hne the 3cr error eUipse for the radio pulsar 47 Tuc U, after 
accounting for the transformation between radio and optical frames and assuming random errors. The optically detected binary companion 
to 47 Tuc U (Uopt) is near the center of the Wll error circle. 
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Fig. 2.— V vsU -V and V vs V - I CMDs for the WF2 chip. The 47 Tuc U MSP companion Uopt is denoted by 'U', a CV by '1' and 
other blue stars are circled. The position of Uopt is slightly redder and brighter than those of the He WDs in NGC 6397 reported in Cool et 
al. (1998). The thicker lines show Bergeron, Wesemael, & Beauchanip (1995) models for 0.5 M0 (dashed line) and 0.6 Mq (solid line) CO 
WDs. The thinner lines are SAROl models, for He WD masses of 0.406, 0.360, 0.327, 0.292, 0.242, 0.196 M©, with mass decreasing towards 
the red, and 0.169 Mq (dot-dashed line). Cooling ages are as shown. These models have been plotted assuming that (m — M)o = 13.27, 
the mean of the 9 distance modulus values reported in Zoccali et al. (2001), E{B — V) = 0.055 (Zoccali ct al. 2001), Au/Ay = 1.51 and 
Aj/Av = 0.60 (Holtzman ct al. 1995). Using (m - M)q = 13.27 the 47 Tuc distance is 4.5 kpc. 
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Fig. 3. — Power spectra and phase plots for Uopt. The upper panel shows power spectra {V- and 7-band, as labeled) for Uopt, with power 
shown in Lomb-Scargle units (Scargle 1982). The binary frequency of 47 Tuc U measured from radio observations is indicated with an arrow. 
The lower panel shows phase plots for the V- and (7-band data for Uopt, where fractional intensities are plotted and the errors are from the 
time series rms values scaled by the square root of the number of points in each bin. 'F' denotes the forced fit and 'U' the unforced fit (see 
§2.3 for description). The point near phase 0.5 in the [/-band data is formally only significant at the 2.3o- level. No evidence is found for 
variation in the /-band, and so the phase plot is not shown here. 



